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It is argued that the experimentally observed baryon stopping indicates a non-monotonous be- 
haviour as a function of the incident energy of colliding nuclei. This can be quantified by a midrapid- 
ity reduced curvature of the net-proton rapidity spectrum and reveals itself as a zigzag irregularity 
in the excitation function of this curvature. The three-fluid dynamic calculations with a hadronic 
equation of state (EoS) fail to reproduce this irregularity. At the same time, the same calculations 
with an EoS involving a first-order phase transition and a crossover one into the quark-gluon phase 
do reproduce this zigzag behaviour, however only qualitatively. 



o 

CN ; I. INTRODUCTION 

a , 

, ^ . A degree of stopping of colliding nuclei is one of the basic characteristics of the collision dynamics, which determines 
a part of the incident energy of colliding nuclei deposited into produced fireball and hence into production of secondary 
particles. The deposited energy in its turn determines the nature (hadronic or quark-gluonic) of the produced fireball 
and thereby its subsequent evolution. Therefore, a proper reproduction of the baryon stopping is of prime importance 
for theoretical understanding of the dynamics of the nuclear collisions. 

A direct measure of the baryon stopping is the net-baryon rapidity distribution. However, since experimental 
information on neutrons is unavailable, we have to rely on proton data. Presently there exist extensive experimental 
Q , data on proton (or net-proton) rapidity spectra at AGS and SPS energies. These data were analyzed 

^ ' within various models (lol-[l8 1 The most extensive analysis has been done in [l4l . Il7j . Since that time new data at 
'. SPS energies have appeared S-Q- Therefore, it is appropriate to repeat this analysis of already extended data set. 
Here it is done within the framework of the model of the three-fluid dynamics (3FD) The 3FD model with the 
hadronic EoS [l^ reasonably reproduces a great body of experimental data in a wide energy range from AGS to SPS, 
^. seelHll^il. 

a^ 

O ■ II. ANALYSIS OF EXPERIMENTAL DATA 

(N , 

' Available data on the proton (at AGS energies) and net-proton (at SPS energies) rapidity distributions from central 
, heavy-ion collisions are confronted to results of the calculations within the 3FD model with hadronic EoS [l^ in right 
^ I ' panel of Fig. [1] Only the midrapidity region is displayed in Fig. [TJ since it is of prime interest in the present 
. . [ consideration. The data at lOA GeV are repeated in the SPS-panels in order to keep the reference spectrum shape 
^ ■ for the comparison. The data are plotted as functions of a "dimensionless" rapidity (y — ycm)/ycrm where ycm is 
k>( [ the center-of-mass rapidity of colliding nuclei. In particular, this is the reason why the experimental distributions 
. are multiplied ycm- This representation is chosen in order to make different distributions of approximately the same 
' width and the same height. This is convenient for comparison of shapes of these distributions. 

As seen from Fig. [U description of the rapidity distributions with the hadronic EoS is reported in [13, [13 ■ The 
reproduction of the distributions is quite good at the AGS energies and at the top SPS energies. At AOA GeV the 
description is still satisfactory. However, at 20A and 30A GeV the hadr.-EoS predictions completely disagree with 
the data [l^- At 20 A GeV instead of a bump at the midrapidity the hadronic scenario predicts a quite pronounced 
dip. In order to quantify this discrepancy, it is useful to fit the data by a simple formula 

dN 

— — = a(cxp{-(l/ti;s)cosh(2/ - j/cm - ys)} + cxp {-(l/iUs) cosh(y - ycm + ys)}) (1) 
ay 

where a, y^ and are parameters of the fit. The form ([l} is a sum of two thermal sources shifted by ±j/s from the 
midrapidity. The width Ws of the sources can be interpreted as Ws = (temperature) /(transverse mass), if we assume 
that collective velocities in the sources have no spread with respect to the source rapidities ij/s. The parameters of 
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FIG. 1: Rapidity spectra of protons (for AGS energies) and net-protons (p — p) (for SPS energies) from central collisions of 
Au+Au (AGS) and Pb+Pb (SPS). Experimental data are from collaborations E802 [H, E877 01, E917 [3|, E866 0], and NA49 
[^-0]. The percentage shows the fraction of the total reaction cross section, corresponding to experimental selection of central 
events. Left panel: Solid lines represent calculations within the 3FD model with hadronic EoS [l9[. Right panel: Solid lines 
connecting points represent the two-source fits by Eq. (HJ. The dashed line is the fit to old data on Pb(158A GeV)-|-Pb [3, 
these data themselves are not displayed. 



the two sources are identical (up to the sign of ys) because we consider only collisions of identical nuclei. Results of 
these fits are demonstrated in the right panel of Fig. [T] 

The above fit has been done by the least-squares method. Data were fitted in the rapidity range \y — ycm\/yc7n < 0.7. 
The choice of this range is dictated by the data. As a rule, the data are available in this rapidity range, sometimes 
the data range is even more narrow (40 A, 80 A GeV and new data at 158 A GeV @). We put the above restriction in 
order to treat different data in approximately the same rapidity range. Notice that the rapidity range should not be 
too wide in order to exclude contribution of cold spectators. 

We met problems with fitting the data at 80 A GeV [1] and the new data at 158 A GeV Q. These data do not 
go beyond the side maxima in the rapidity distributions. The fit within such a narrow region results in the source 
rapidities y^ very close (at 80 A GeV) or even exceeding (at 158 A GeV) ycm and a huge width Ws- As a result, the 
normalization of the net-proton rapidity distributions, as calculated with fit ([1}, turns out to be 330 (at 80 A GeV) 
and 400 (at 158A GeV), which are considerably larger than the total proton number in colliding nuclei (=164). To 
avoid this problem, we performed a biased fit of these data. An additional condition restricted the total normalization 
of distribution ([T]) to be less than the total proton number in colliding nuclei (=164). This biased fit is the reason 
why the curve fitted to the new data at 158A GeV does notperfectly hit the experimental points. In particular, 
because of this problem we keep the old data at 158A GeV Q in the analysis. We also use old data at AOA GeV, 
corresponding to centrality 7% [1], instead of recently published new data at higher (5%) centrality since the data 
at the neighboring energies of 20A, 30A and 80A GeV are known only at centrality 7% @. Similarity of conditions, 
at which the data were taken, prevents excitation functions, which are of prime interest here, from revealing artificial 
irregularities. 



III. PARAMETERS OF THE FIT 



Energy dependence of parameters ys and Ws deduced from these fits of the data revels no significant irregularities: 
they monotonously rise with the energy. At the same time, inspection of the evolution of the spectrum shape with 
the incident energy rise reveals an irregularity. Beginning from the lowest AGS energy to the top one the shape of 
the spectrum evolves from convex to slightly concave at lOA GeV. However, at 20A GeV the shape again becomes 
distinctly convex. With the further energy rise the shape again transforms from the convex form to a highly concave 
one. In order to quantify this trend, we introduce a reduced curvature of the spectrum in the midrapidity defined as 
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FIG. 2: Midrapidity reduced curvature (left panel) and midrapidity value (right panel) of the (net)proton rapidity spectrum as 
a function of the center-of-mass energy of colliding nuclei as deduced from experimental data and predicted by 3FD calculations 
with hadronic EoS (hadr. EoS) as well as with EoS's involving a first-order phase transition (2-ph. EoS) and a crossover 
phase transition (crossover EoS) into the quark-gluon phase [2^. The thin long-dashed line on the right corresponds to the 
hadr. -EoS calculation without fragment production, i.e. without coalescence. 



follows 

/ o (fN\ , ( dN\ , , ^2 ^ . , 2 , N 

(^y = I I / [Vcm— I = KVcm/Ws) [smh ys -WsCOshys) . (2) 

V y / y=ycm \ ^ y=ycm 

This curvature is dcfirLcd with respect to the "dimensionless" rapidity (y — ycm)/ycm- The factor 1/ {ycmdN/dy)y^y 
is introduced in order to get rid of overall normalization of the spectrum, i.e. of the a parameter in terms of fit ([1]). 
The second part of Eq. ([2]) presents this curvature in terms of parameters of fit ([1]) . 

Thus, the reduced curvature, Cy, and the midrapidity value, {ycmdN/dy)^^^ , are two independent quantities 
quantifying the the spectrum in the midrapidity range. Excitation functions of these quantities Values of the curvature 
Cy deduced both from fit ([T|) to experimental data and from the same fit of the results of the 3FD calculation with 
different EoS's are displayed in Fig. [2] 3FD calculations with hadronic EoS (hadr. EoS) [l^, as well as with EoS's 
involving a first-order phase transition (2-ph. EoS) and a crossover phase transition (crossover EoS) into the quark- 
gluon phase are presented. Notice that a maximum in ycm{dN/dy)cm at s^/^ = 4.7 GeV happens only because 
the light fragment production becomes negligible above this energy. The 3FD calculation without coalescence (i.e. 
without the fragment production) reveals a monotonous decrease of ycm{dN/dy)cm beginning from 5^/2 ^ 2.7 GeV, 
i.e. from the lowest energy considered here. 

To evaluate errors of Cy values deduced from data, we estimated the errors produced by the least-squares method, 
as well as performed fits in different the rapidity ranges: \y — ycm\/ycm < 0.5 and \y — Vcml/ycm < 0.9, where it is 
appropriate, and also fits of the data at 80^ GeV Q and the new data at 158^ GeV Q with different bias on the 
overall normalization of the distributions: A^prot. £ 208 (i.e., half of the net-nucleons can be participant protons) 
and A'^piot. < 128 (which is the hydrodynamic normalization of the distribution). The error bars present largest 
uncertainties among mentioned above. The lower point at s^^^ = 17.3 GeV corresponds to the new data at 158A 
GeV. Its upper error, as well as that of 80A-GeV point, results from the uncertainty of the normalization. 

The irregularity observed in Fig. [l]is distinctly seen here as a zigzag irregularity in the energy dependence of Cy. 
A remarkable observation is that the Cy curvature energy dependence in the first-order-transition scenario manifests 
qualitatively the same zigzag irregularity (left panel of Fig. [2]), as that in the data fit, while the hadronic scenario 
produces purely monotonous behaviour. This zigzag irregularity of the first-order-transition scenario is also reflected 
in the midrapidity values of the (nct)proton rapidity spectrum (right panel of Fig. [5]). As for the experimental data, it 
is still difficult to judge if the zigzag anomaly in the midrapidity values is statistically significant. In the conventional 



representation of the data without multiplying by j/cm, the irregularity of the {dN/dy)cm data is hardly visible [24 
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FIG. 3: Dynamical trajectories of the matter in the central box of the colliding nuclei (4fmx4fmX7cm4fm), where 7cm is the 
Lorentz factor associated with the initial nuclear motion in the cm. frame, for central {b — 0) collisions of Au+Au at 
and lOA GeV energies and Pb+Pb at 20A GeV. The trajectories are plotted in terms of baryon density (ub) and the energy 
density minus ub multiplied by the nucleon mass (e — niNnB)- Only expansion stages of the evolution are displayed for two 
EoS's. Symbols on the trajectories indicate the time rate of the evolution: time span between marks is 1 fm/c. 



The crossover EoS represents a smooth phase transition, therefore, it is not surprising that it produces only a wiggle 
in Cj,, which again only qualitatively resembles the experimental behavior. 

The baryon stopping depends on a character of interactions (e.g., cross sections) of the matter constituents. If 
during the interpenetration stage of colliding nuclei a phase transformation^ of the liadronic matter into quark- 
gluonic one happens, one can expect a change of the stopping power of the matter at this time span. This is a 
natural consequence of a change of the constituent content of the matter because hadron-hadron cross sections differ 
from quark-quark, quark-gluon, etc. ones. This naturally results in a non-monotonous behaviour of the shape of the 
(net)proton rapidity-spectrum at an incident energy, where onset of the phase transition occurs. 

However, if even the same friction is used in both phases, the calculated (with 2-ph. EoS) reduced curvature still 



reveals a zigzag behaviour but with considerably smaller amplitude, as it was demonstrated in [2J]. This happens 
because the EoS in a generalized sense of this term, i.e. viewed as a partition of the total energy between kinetic and 
potential parts, also affects the stopping power. The friction is proportional to the relative velocity of the counter- 
streaming nuclei [iTj . Therefore, it is more efficient when the kinetic-energy part of the total energy is higher, i.e. 
when the EoS is softer. This is precisely what the phase transition does: it makes the EoS essentially softer in the 
mixed-phase region. The latter naturally results in a non-monotonous evolution of the proton rapidity spectra with 
the energy rise. 

Figure [3] demonstrates that the onset of the phase transition in the calculations indeed happens at top- AGS-low- 
SPS energies, where the zigzag irregularity takes place. Similarly to that it has been done in [28|, the figure displays 
dynamical trajectories of the matter in the central box placed around the origin r = (0,0,0) in the frame of equal 
velocities of colliding nuclei: |x| < 2 fm, \y\ < 2 fm and |z| < 7cm 2 fm, where 7cm is Lorentz factor associated with 
the initial nuclear motion in the cm. frame. Initially, the colliding nuclei are placed symmetrically with respect to 
the origin r = (0,0,0), z is the direction of the beam. The e-Ub representation is chosen because these densities 
are dynamical quantities and, therefore, are suitable to compare calculations with different EoS's. Subtraction of the 
misinB term is taken for the sake of suitable representation of the plot. Only expansion stages of the evolution are 
displayed, where the matter in the box is already thermally equilibrated. The size of the box was chosen to be large 
enough that the amount of matter in it can be representative to conclude on the onset of the phase transition and 



^ The term "phase transition" is deliberately avoided, since it usually implies thermal equilibrium. 
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to be small enough to consider the matter in it as a homogeneous medium. Nevertheless, the matter in the box still 
amounts to a minor part of the total matter of colliding nuclei. Therefore, only the minor part of the total matter 
undergoes the phase transition at lOA GeV energy. As seen, the trajectories for two different EoS's arc very similar 
at AGS energies and start to differ at SPS energies because of the effect of the phase transition. 



In conclusion, it is argued that the experimentally observed baryon stopping indicates a non-monotonous behaviour 
as a function of the incident energy of colliding nuclei. This reveals itself in a zigzag irregularity in the excitation 
function of a midrapidity reduced curvature of the (nct)proton rapidity spectrum. The energy location of this anomaly 
coincides with the previously observed anomalies for other hadron-production properties at the low SPS energies 
j25l [29} . The 3FD calculation with the hadronic EoS fails to reproduce this irregularity. At the same time, the same 
calculations with the EoS involving a first-order phase transition (within the Gibbs construction) and a crossover one 
into the quark-gluon phase (23| reproduce this zigzag behaviour, however only qualitatively. Preliminary simulations 
with the EoS of [sO] , also based on the first-order phase transition but within the Maxwell construction, show the same 
qualitative trend. It is argued that the non-monotonous behaviour of the baryon stopping is a natural consequence of 
a phase transition. The question why these calculations do not qualitatively reproduce the zigzag irregularity deserves 
special discussion elsewhere. It is very probable that the Gibbs and Maxwell constructions are inappropriate for the 
fast dynamics of the heavy-ion collisions [3l|, HH and metastable states should be explicitly considered. 



It is somewhat suspicious that the zigzag irregularity happens at the border between the AGS and SPS energies. It 
could imply that this irregularity results from different ways of selecting central events in AGS and SPS experiments. 
Moreover, all data at SPS, except for those at the energy 158A GeV, still have preliminary status. Rapidity range 
of data at 80A and 158A GeV is too narrow for their reliable analysis. All this indicates that upgraded data are 
desperately needed to either confirm or discard the existence of the zigzag irregularity. It would be highly desirable 
if these new data are taken within the same experimental setup and at the same centrality selection. Hopefully such 
data will come from new accelerators FAIR at GSI and NICA at Dubna, as well as from the low-energy-scan program 
at RHIC. In this respect the NICA project looks the most promising, since it covers the whole energy range of interest. 
The discussed possibility to detect neutrons in the NICA/MPD detector could provide information on real baryon 
stopping irrespective of any assumptions on redistribution of baryon charge between net protons and neutrons. 
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